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NUMERICAL EVALUATION OF THE WAKE-SURVEY EQUATIONS 
FOR SUBSONIC FLOW INCLUDING THE EFFECT 
OF ENERGY ADDITION 


By Donald D. Baals and Mary J. Mourhess 
SUMMARY 


Direct-reading tables and charts are presented for 
determining the drag or thrust coefficients from wake- 
Survey measurements in the subsonic speed range. For 
flows wherein no energy is added, the point drag coef- 
PECHeML LS ShOwn bO) be) an explicit function, of ithe stream 
Mach number Mp), the static-pressure coefficient at the 
wake station P;, and the total-pressure-loss coeffi- 
cient AH/a55 where AH is the total-pressure loss 


and do is the stream dynamic pressure. Values of the 


point drag coefficient are tabulated for a wide range 

of values of these parameters. Inasmuch as the tabulated 
coefficients (either drag or thrust) represent the point 
values, Wiich are independent of the integration of the 
wake, the charts or tables in the form presented are 
general in application. 


For flows wherein energy is added, such as flows 
behind propellers or heated radiators, an additional 
parameter, which is a function of the stagnation- 
temperature rise, must be considered. Values of the 
point dras coeifiecient that include the effects of the 
addition of energy are tabulated. 


INTRODUCTION 


ar 


In the field of aeronautical research, wake pressure 
surveys have been used increasingly for determining 
profile drag, internal drag, jet thrust, and related 
factors. The evaluation of the drag or thrust from 
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pressure surveys is based on the solution of the momentum 
equation. Tor hish=sneed flows, in which the are musi be 
considered compressibie, the evaluation of the drag or 
thrust must include the variation in the density of the 
air; this variation in the density necessarily complicates 
the solution of the momentum equation. For flows wherein 
energy has been added, the density variation is of 
importance even for the low-speed conditions. 


Because ‘a Maree: number or pressure) readineas) ane 
Usvaliv weculced ho make eaicomplete Wakes Suleiveiy aimee 
essential that solutions of the wake-survey equations 
iilone Bill Taos) Ope “2c ILOMy lols jose Sioiagelel alial) 2) “aE Toga) THe) joe C 
rapid and accurate computation. Various simplifying 
techniques have been Geveloped for evaluation of the 
drags coefficient from wake pressure measurements by 
means) of charts or Gables, bub these methods Nave requames 
EXCESS Vel internolatvon OF Compubauaom Cue) To: tie yor 
on the paraneters employed. “When Simpl Tt hyan ea asisump Eases 
as to the shave of the wake profile ana the constancy of 
the Statice pressure across the wake are made, intesracon 
techniques may be combined with the solution of the wake-= 
survey equations (references 1 and 2) to reduce greatly 
the recurred compuuation. | This) meunodliias) been a eune 
uUseiuie im Ceverminine The ‘Sec lon prods 6 cia IC OS iiememomine 
OL alriouds however, such a methodmaist now general) am 
application and is limited by the opiginal assumptions 
of uniform static pressure and a definite wake profile, 


The method found to be most generally applicable at 
the Langley 8-foot high-speed tunnel consists in the 
point-by-point evaluation of the drag coefficient by 
means of direct-reading charts with the final integra- 
tion jerhormed in a weparace, Step acd liye m lenraiete 
li nal oval Moje Piclalen (cligkeyer @ ele ibsP weal Sratiny VaLis) Moo allel’ tiny Sulake” wwlale 
point-by-point values can be determined from simple 
charts, The charts require no assumptions as to wake 
jOleonial INS) ona Monaesisibhaien eeaVlaleial missy | Pays) ac llieecjoibewieny eSCT c's OC. 
for application to both two- and three-dimensional flows 
1s provided by performing the integration in a Sseparave 
step. In this report the solution of the wake equations 
has been developed for both isoenergic flows (constant 
total energy) and flows wherein energy has been added; 
thus the equations are generally applicable to flows 
behind, airfoils, propellers, and radiators, . The funda= 
mental principles may be applied to the determination of 
the thrust from jet units; however, consideration must be 
given to the change in the value of the ratio of specific 
heats and the momentum of the fuel. 
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SYMBOTS 
a speed of sound, feet per second 
A area, square feet 
Ap reference area, square feet 
D 
Ch drag coefficient 
IoAp 
c chord, feet 
Cg section profile-drag coefficient 
Cg! point drag coefficient 
Cc specific heat at constant pressure (for air 
Pp f 3 
0.2) Btu/1lb/oF) 
D drag, pounds 
E energy added, foot-pounds per second 
Fi compressibility factor oe ) 
q 
g acceleration of gravity (32.2 ft/sec@) 
H tOval pressure, pounds per square foot 
AH Ho = 
er total-pressure-loss coefficient [SS 
J mechanical equivalent of heat (778 ft-1b/Btu) 
E 
K energy-input factor oS 
Je,,gmly 

M Mach number (2) 
m mass flow rate, slugs per second (pAV) 
p Static pressure, pounds per square foot absolute 

Paes 
oy static-pressure coefficient wre RUT 

©) 
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q dynamic pressure (tov?) 

T static Pe ee OF absolute 

ay stagnation temperature, CF absolute 

aM stagnation-temperature rise, °F oe = Ta 
V VeLOCiiy, (keeu pee Second 

y distance across wake, feet 

0 Gensity, Slugs per cubse foou 

¥ ratio of specific heats (for air, 1.100) 


Subscripts 6, 1, and 2 reler vores son Stacroms Cec uae 
Maceq am) weemiee le 


DISCUSSION OF THHORY AND METHODS 


Fundamental Relations 


The detailed solutions of the basic wake-survey 
equations for flows with and without the addition of 
energy have been presented in reference 3. A theoretical 
analysis) Of (ENG Dasie ASSuUMpinL ONS) requarcd) Tory ties ole 
tion of the wake-survey equations for isoenergic flows 
has been presented in reference }. 


The basic form of the drag equation as derived from 
the momentum relation is 


n 


D = | p1V,(V > va) oa, 
wake 
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oe Py a, ‘ 
oe | sre \We.7 Ve) aay 
oO" - 
wake 


ce of MS TN Peo Mitek, Mey cles, 

b AR 46 ar Pe Po Po Vo 
wake 

3 / ehxeom'G (2) 

AR *) qo 
wake 


With the usual assumption that the total pressure at sta- 
tion 2 (where po = Da! is equal to the total pressure 

at station 1, the drag coefficient can be evaluated as 
indicated in the appendix. The numerical solution of 
equation (1), however, is difficult and requires excessive 
computation for direct use. 


-_— 


1) 


For convenience in presenting and discussing the 
solution of equation (1) the point drag coefficient is 
defined as 


“VOW OO) * 


Negative values of the point drag coefficient indicate 
that the system is producing thrust. 
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Isoenergic Flows 


Isoenergic flows are found over aerodynamic bodies 
or through internal-flow systems when no avpreciabls 
chenge in stagnation temperature occurs. An analysis 
of the terms constituting the point. drag) cocificient sor 
isoenergic flows indicates that the point drag coeffi- 
cient is an explicit fumction of the free-stream Mach 
number M., the static-pressure coefficient in the 
wake Pj, and the totval-pressure-loss coefficient AH/q,- 


The pressure coefficients wsea are merely an expression 
or the measured Gata in coetiictent worm. Because sac 
values of these parameters can be easily determined from 
the test data, their use provides a convenient and direct 
method for evaluation ef the point crag coerficient. As 
CzxplLiGlEyexpression Rom Lhe pornit, Geas. Coearietent, am 


terms of the parameters M,, P,, and AH/q, is com- 


plicated; but the coefficient can be easily cetermined 
by usuai methods for given values of these parameters 
(see appendix). 


Values of the point drag coefficient for a wide 
range of values’ of pressure coefficients’ P, and A/a, 
and for given values of stream Mach number M, are 
presented in table I. The range of total-pressure-loss 
coefficients has been extended into the negative region 
(which indicates a total-pressure increase) to permit 
evaluation under conditions of low energy input approaching 
isoenergic flow. The application of table I to flows 
wherein energy has been added will be discussed in a 
later section. A large range of positive values of Pz 
has been included to permit computation of internal duct 
flows for which the static-pressure coefficient may 
approach unity. 


Other forms of the parameters and other methods of 
presentation than those used in table I may be more 
suitable for certain applications. Instead of Mach 

Eg 


Ho = Po Ho - Po 
3 


number M the pressure ratios 


“ 


oO? = 
Po ae Ho 
or “ may be used; the pressure coefficients Ps 
ie) r 
= = H Ta H 
AH : Pa Se fe) 1 
and =— ‘may be replaced by = and Sail 
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by = DP Hea; 
ub Off ats fe) ne 


Ho - Po Hy - Py 
have an advantage when they can be determined directly 
from the measured wake and stream pressures. Tables and 
Charts.of the point drag coefficient for isoenergic flow 
in terms of these two parameters are presented in refer= 
ence 5 with charts and tables for rapid evaluation by an 
approximate method. 


respectively. The two parameters 


Flows Wherein Energy Is Added 


The fundamental relations expressed in equation (1) 
are. Correct Tor the evaluation of the drag or’ thrust 
Goctficient. for flows wherein energy has’ been added, such 
as flows through radiators or propellers. ‘Tre evaluation 
of the density ratio Po/Pos however, involves an added 


parameter that is a function of the energy input. The 
evaluation of the density ratio P5/P, as developed in 


appendix B of reference 4 becomes 


Mies or 22 
1 + M oc 
p 
re = : ce) 
Po oe a E 
ne Moo + 
2 Je, gt, 


tie erder to determine the dras cosiiicient from equa= 
tion (1) for flow conditions wherein energy has been 


added, an additional parameter eee ah must be deter- 
Jc,gmT, 


mined; this parameter is deadematen | k 


An inspection of the terms constituting the point 
drag coefficient oie )) indicates that the 


p 
gape  ) @ can be expressed as a function 


oF ii. z —.,. The last part of equation (2) 
fo) 
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aay =— Gan be expressed as a function of Ilys 


—,, , ang, K. (The two, quantities of which they prodie as 


the point drag coefficient can therefore be determined 
independently. Values of the two parts of equation (2) 
required to evaluate the point drag coefficient for flows 
wherein energy has been added are presented in tables If 
and III. 


The value of the energy parameter K can be com= 
puted directly for flow. conditions where the energy, 
input and mass flow are measurabie’ and are uniform across 
the survey plane. Such a condition might be attained 
for fiow through an efricient radiator installation. | Hor 
flow conditions as found behind a propeller, however, 
where the energy input and elemental mass fliow are not 
uniform and are difficult to measure, an evaluaticn of 
the energy parameter can better be made by experimental 
methods. 


The energy equation from the free stream to the wake 
station becomes 


2 


Oo 
Seg TEL SM se 
2 peo 


For stagnation conditions at the free-stream and wake 
stations, the energy parameter becomes 


Jagei ts 


Hoke ! 
tT Ty 


fo) 
AVE? 


—— 
=) 


where the prime refers to the stagnation temperature, 
The energy parameter becomes, then, merely the ratio of 
the stagnaticon-temperature rise to the absolute stream 
static temperature. 
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Inasmuch as the stagnation-temperature differ- 
ence 1)! - TT)! is small and must be measured accurately, 


Stagnation-temperature thermocouples or resistance 
thermometers may be connected and calibrated to read 

the stagnation-temperature difference directly. The 
free-stream stagnation-temperature reference may be 
obtained by installing one of the temperature-measuring 
elements outside the wake. The determination of the 
absolute free-stream static temperature does not require 
sreat accuracy; therefore, conventional methods may be 
used. 


USE OF TABLES AND CHARTS 


The values presented in tables I to III for evalua- 
tion of the point drag coefficient are difficult to apply 
directly because interpolation is required. The tabulated 
values have therefore been plotted in the form of direct- 
reading charts. 


Tsoenergic Flows 


Examples of the types of plot found useful for 
Evaluation of the point drag coefficient for isoenergic 
flow are presented in figure 2 for M, = 0.20 and 0.70, 
respectively. A representative cross plot of the tabular 
data of table I is presented in figure 3. It should be 
pointed out that the scale of the plots of point drag 
coefficient presented in this report is too small for 
accurate work. These plots have been included, however, 
to indicate the general form of the curves and to provide 
a basis for more accurate large-scale plots based on the 
original tabular data. 


The charts in their present form are especially 
applicable for use in high-speed wind tunnels, where 
tests are usually run at fixed values of stream Mach 
number, and the value of the static-pressure and 
total-pressure-loss coefficients can be determined 
directly from the pressure records and wind-tunnel 
calibration. For use at values of stream Mach numbers 
different from those tabulated, the chart for the 
nearest tabulated Mach number may be used with but small 


1 : 
error +15 percent jie For sreater Bsecursey “abt (the 3nter= 


mediate Mach numbers a linsar interpolation may be 
assumed (fig. 3). 
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The presentation of the tabulated results may be 
Changed ifromethat used in Figure 2°te fit .various con= 
Cq! 

AH 
/2., 


Gitvons. Om exeanyve. a6 le nie auae may be plotted 
instead of cg! in order to attain greater accuracy on 

a small-scale plot. For many tests the wake static- 
pressure coerficient remains essentially constant with 
changes in Mp; therefore the point drag coefficient 

may be plotted against AH/q for a range of values 

» the use of the tabutaged 


data in this form may be applicable to flight tests, 


fe 


of M, at given values of P 
1 


Flows Wherein Energy Is Added 


In figure h are piteSsenced —represcntabi ve jplous sor 
tables Ifvand Iii for evaluating, the, point drag cogkir= 
Gient, for Lows,where in ieneseyeus addedsy « dln iss) oie pues 
shows that an increase in stagnation temperature of the 
order of 19 F (K = 0.002) can produce a significant 
effect on the evaluation of the drags coefficient at Low 
values of AH/Q5- For many cases, however, the energy 


effect may be considered negligible = as; for example, 
the flow behind a lightly loaded propeller - and the 
tabulated results for isoenergic flow (table I) may be 
used directly. A comparison of the numerical result for 
Ghe Sexperimentals yale Oi 9 twas “cha ome — a0 
indicates the magnitude of the energy effect and which 
form ‘ofthe’ equations “should be “used sy" The” iscenersare 
values of cq! presented in table I correspond to the 


Values of © cg" ror kK ="0 “presented sn tables fi aad ie 


Because the solutions of the equations for usoenersic 
flow Can. be /GeLrernined - trom Gale) Podirecily aa one 
step, table I should be used whenever the energy effect 
can ~oe, mesieeted. 


It shoulc be pointed out that for flows wherein 
energy is added, certain combinations of the parameters 


Mos. Pale = ang) K) specify ‘supersonic flow in the 
fe) 


wake. For such flows the measurement of the wake pres- 
sures is subject to corrections that are dependent on 

the intensity of the shock at the local supersonic speed. 
No attempt has been made in this report to evaluate 

these shock corrections to the measured pressures. 
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Integration Techniques 


The evaluation of the total-drag coefficient involves 
the integration of the point-drag-coefficient profile. 
Tnasmuch as the evaluation of cau dus| PUalslejclevalalevalin (Os 
the intesration process, great flexibility in the appli- 
Cation of the proposed wake-survey techniques to various 
types of flow can be attained. The flows behind air- 
foils, within internal-flow systems, or through pro- 
pellers may therefore be calculated by similar methods 
that vary only in the integration technique employed. 
The wake profile and the manner in which the wake is 
surveyed will determine the optimum type of integration. 


With the values of the point drag coefficient deter- 
maned from the charts or tables presented in this report, 
the value of the total-drag coefficient becomes 


aihorli 
R 
wake 


where A, is the reference area upon which Cp is to 


be based. The integration technique employed should be 
determined from considerations of expediency and the 
accuracy of integration desired. 


For two-dimensional flows behind airfoils the 
airfoil-section drag coefficient becomes 


where the integral represents the area under the cg!- 
curve when expressed in the same units as the chord c. 
NUMERICAL SXAMPLES 
Isoenergic Flow 
Assume the following set of conditions, which might 


correspond to the center of the wake of an airfoil at 
high speeds: 
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M, = 0.70 
EL = Oe 
Au = One 
do 


An interpolation of table I or the use of the large- 
scale plots of table I (fig. 2(b), for example) gives 
the value of the point drag coefficient ecg! of 0.196. 
The resulting positive value indicates drag. 


Flows Wherein Energy Is Added 
(1) Assume the following set of conditions, which 


might correspond to the flow behind a highly loaded pro-= 
peiiler wom wae Gilad CGCiacliipaenm: 


My = 0,20 

Er = G25 

AH 

SB E26) 56 
do 


(The negative value of AH/4, indicates total-pressure 
increase. ) 


Ant = le Oo 


T= 500° F aps. 


From the value of absolute static temperature and 
stagnation-temperature rise, 


I,.0 
500 


K = 


0.008 


AM Ante rpoLabion of tables TE and Lif or the use of the 
large-scale plots of tables II and III (fig. 4, for example) 


WACA ARR No. L5H27 13 


gives the following values of the parts of equation (2) 
Ghat constitute the point drag coefficient: 


ae 1/2 
) a ne 2.2368 
Po ae 


(2) I ED ass 
ce ce 


The value of the point drag coeificient thus becomes 


Ga (2.2760) (= Ora217) 


~0.4959 


Hi} 


whe resulting negative value indicates thrust. 


(2) Assume the following set of conditions, which 
might correspond to the flow behind a radiator measured 
near the duct outlets: 


H== 
iF 
= 600 Btu/sec 
gem = 30 lb/sec 
eS gO° raps’, 
We 1007 
AH 
== Oat 
sro 


ian = Oem 
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Now 


J ChE gullies 
600 
O.2h x 50 x 150 


i] 
io) 
kh 

\O 
t= 


From these values of the parameters, 


py\\/2/¢ ily 2 
ayia = 199.09 


and 
pen 72 qaniy2 
(=) - S = 0.0889 
Po/ 8g 


The value of the point drag coefficient thus becomes 


= 1.5503 x 0.0889 


OmtS73 


The resulting positive value indicates drag. 
Langley Memorial Aeronautical Laboratory 


National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 


EQUATIONS FOR COMPUTING TABULATED VALUES 


OF DRAG COEFFICIENT 


The various terms of the drag-coefficient edquation 
(equation (1))can be expressed by the following relations. 


MEecmew! a he polnted out “sheatt "oy dering tion ba = De 


and by assumption H, = Ho. (A Getailed derivation of 
the following terms can be found in appendix B of refer- 
ence 3.) 


/ 2 
=) 2 7 i Fg \L/e 
HW = 75 7 
do HO a ) 


Ne | as W/Z) iL 
a\ fi By yy Fo, /2 
@) 


Py re) ay 


1 + Mo? -—= 

Ss i ae 
es) 21 

% a Mo- + K 


Where K. is the energy-input factor. (For isoenergic 
pow kK = .0.) 


for the use), of the equations in the foregoing form, 
Some initial value of ene cf the pressures musu be 
assumed. Because the values of the various parts of 
equation (2) are determined fundamentally by the pressure 
Fatio, the numerical value of a pressure has no signifi= 
cance; therefore any positive value may be used. 
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The détailed computing procedure is as follows: 


(1) Given M, and a range of values of Py 
end AH/q + (The value of K must also be considered 
for flows wherein energy is added.) 


(2) In order to determine the free-stream pressures 
corresponding to M,, assume any convenient value of 
stream static pressure p. (1000 lb/sq ft is assumed 
here) and use the following relations: 


a iaelo! 
IB VEE IO) ae Fe do 


where vy = 1,400 and F, is obtained from table IV 
fo) 

fon che, value on Wiis) iiehvalvcy of ain cana som ic 

Gleneagle, wore Swed Waslwes! ow i ebiel 9. ae eu the 

following equation: : 

(ee) 


Me = eae B ie earl 
Np oa Ae NG) 


(4) The step-by-sten computing procedure is as 
follows: 
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Relation 

Sale) 
C7) Fg, from a plot of table IV for value of (6) 
(8) He) iP a= (h) - Do 


ee eas wee) 
Pa Po 


Leb hota Ole Vitor ve auenie Toten (S) 


ee =i 


eee —~ 
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a as Rie Gi a LL Ic a 
Column Relation | 


| 
SFL L fev Gag oit | 
Le) Fics qeae 


y -l G 
1+ ye —S 
| 2 Di OW Ge 
(16) Pathe a 
y 
Mo" + K 


Le Oat (lL) 


@) 
fe 
++ 


-—_— 
fH 
—~J 
Se O 
ay 
4 
SS 
n> 
1 
bt 
ON 
—— 
ae 
po 


ih es 0.2M,* i 


Hoe asiosine cee Tle, = 0 


i) m= es) 


(19) a0 = eS o ay" “| 


awe -) Le 


| 


I 
ne) 
x 
YY 
Sa 
ee 
= 
x 
Se 
CO 
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(L.) The numerical solutions of the various steps of 


the computing equationsare Biven for the following 
examples: 


——S ewe ee 


1000 


i | | 
3143 | 28 33 
a 1487.00, - | 1025.28 1387.08 
| 

i) 102.9 =i ae 
(2) ibaa | UG 
(3) Monae 6. ae * 1007 104), .3 
(L,.) 128),.18), 1042.28 129.88), 
(5) 2h4 02h 35.28 eae 
(6) ee 0.0350 0.208), 
(7) 1.079 i Olain be Oasis) 
(3) 28h .18h 2.28 2h.9.8dh 
(9) 0.282 | 0.0.23 0.2h.99 | 
(10) 1.0960 | 1.0150 1.0815), 
am) 0.7560 | Ila siay OR Sy 
(12) 0.8695 | ey aalee 
(13) aeases Meals 0.5867 
(14) 0.8101 LSC 0.7659 
(15) 1.0145 1.0025 ion 
(16) 0.9762 0.9960 0.82h9 
ie 0.9890 0.9980 0.9083 
(18) 0.1195 | -0.2217 ey 
(19) 011964 | -0:1:959 01137 


Any arbitrary value of py, may be assumed without change 
MirpEic mi nae resuilGe Mom) Enese examples ge = 1000 has 
been used. 


<_e 


Noa 
e 
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TABLE I 


POINT DRAG COEFFICIENT FOR ISOENERGIC FLOW 


fs =o 


8) .6396 


eots 
2s 


O00 7 0 0 D fe) 
° ee ° 


hou OW oo W th tt te) 
. eee ee eo 
iodo i i Dp ce) 


Ooo Too Oh DY Do 0 ©) 
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TABLE I = Continued 
POINT DRAG COEFFICIENT FOR ISOENERGIC FLOW - Continued 
[Mo = 0405] 
Py 
iA =0250] -0.45}] -0.40} -0.35| -0.30] -0.25| -0.20] -0.15| -0.10] -0.05 
% 
0.800} 0.923)| 0.8900/0.8552 0.78 5 ; 6 
2700 “30 ae “yeas yee ig 0 
. poe : iby 5 Fle 354} 6 ae 
‘hoo lg crevalue 98 
300 Bon, 3498 2422 
2250] 22990] .2930| .2868 
2200 “2404 22358| .2310 
189) 2890) 208) 2195 
3075 -0918 50896] .0878 
2050; .0608] .0598| .0588 
2025} 0306} .0300] .0296 
0) ) 0 ) 
=208e 570706 meron = 20298 
=.075 =e eoesaullos gh 
-.100)-.123)|-.121)}-.119) 
-.150/-.1860] -.1830|-.1802 
=.200]-.2)8h | -.2hf6|-.2h,10 
= .250|-.3120|-.307 |-.30 
ee ae oe 
=2500 Bees Boo las 
-.600}- ee = ieee = Ste 
-.700 
-.800 5.3339 1. dBse ae oLot 


0.28),8 
22926 
32664, 


0.2328]0.1650|------|------| ------ |------ | ------ |------] ------ |--~--- |------ 
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